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ABSTRACT 

We investigate the generation mechanisms of MHD Poynting flux in the magnetised solar photo- 
sphere. Using radiative MHD modelling of the solar photosphere with initial magnetic configurations 
that differ in their field strength and geometry, we show the presence of two different mechanisms for 
MHD Poynting flux generation in simulations of solar photospheric magneto-convection. The weaker 
mechanism is connected to vertical transport of weak horizontal magnetic fields in the convectively 
stable layers of the upper photosphere, while the stronger is the production of Poynting flux in strongly 
magnetised intergranular lanes experiencing horizontal vortex motions. These mechanisms may be 
responsible for the energy transport from the solar convection zone to the higher layers of the solar 
atmosphere. 

Subject headings: Sun: Photosphere — Sun: Surface magnetism — Plasmas — Magnetohydrodynam- 
ics (MHD) 



1. INTRODUCTION 

The magnetic coupling between the solar interior and 
the outer layers of the solar atmosphere plays a key 
role in the energy transport between these regions . We 
have recently demonstrated (Shel yag et al. 1 12011 ah that 
vorticity generation in the intergranular lanes of mag- 
netic solar plage region simulations is much more effi- 
cient compared to the non-magnetic case. These vor- 
tices have been detected in observ ations of the pho- 
tosphere and chromosphere (see e.g. iBonet et al.l 2008; 
IWedemever-Bohm fc Rouppe van der Voortll2009[ ). The 
magnetic vortex motions are highly significant as they 
can de posit large amou nts of energy in the solar chromo- 
sphere (|Moll et al.ll2012l ). contribute to the production of 
Poynti ng flux and generate various types of M HD wave 
modes ([Shelyag et al.ll201lU IFedun et al.ll2Qllh . On the 
other hand, using a simulation setup that involved in- 
flux of weak horizontal magnetic field thr ough the bot- 
tom b oundary of the simulation domain, Ste iner et al.l 
(2008) have shown that the Poynting flux is generated 
through the advection of horizontal magnetic fields by 
vertical upflows in the up per, convectively s table layers 
of the solar photosphere. lAbbett fc Fisherl (j2012h have 
recently demonstrated the presence of Poynting flux in 
their radiative MHD simulations of an open-flux region 
of the solar atmosphere. 

In this paper we use forward modelling techniques to 
examine the generation of Poynting flux in a model of the 
solar magneto-convection with different initial magnetic 
field strengths and geometries. Our model includes the 
wavelength-dependent radiative transport and the par- 
tial ionisation of a plasma with solar chemical compo- 
sition. We perform simulation runs with initial verti- 
cal magnetic fields of 10 G and 200 G, as well as with 
changing polarity with unsigned flux of 200 G and with 
weak horizontal magnetic field of 10 G, advected into 
the numerical domain by convective upflows. The effec- 
tiveness of Poynting flux generation is analysed in these 
simulations and we find that Poynting flux generation by 



vortices, outlined by Shel yag et al.l ([2011bt h is more pro- 
nounced irrespective of the magnetic configuration. We 
use a mixed-polarity region to study the dependence of 
Poynting flux generation relative to the degradation of 
activity. 

2. SIMULATION SETUP 

Our numerical modelling of the solar magnetised pho- 
tosphere is undertaken using the MURaM radiative 
MHD code (IVogler et aTTl2QQ5h . The code is described 
in detail in IVoglerl ((2003) and it has been successfully 
used for a wide variety of applications in solar physics. 
Here we provide only a brief description of the code. It 
uses 4-th order central difference operators to calculate 
spatial derivatives, a 4-th order Runge-Kutta scheme to 
advance the solution in time, and hyper diffusive oper- 
ators to stabilise the solution against numerical insta- 
bilities. A non-grey radiative source term is included in 
the MHD equations to account for wavelength-dependent 
radiative energy exchange in the simulated photosphere, 
and a non-ideal equation of state is used, which takes 
into account partial ionisation of the 11 most abundant 
elements in the solar photosphere. Side boundary con- 
ditions are periodic, and the top boundary is closed for 
in- and outflows, while allowing horizontal motions of 
fluid and magnetic field lines. The bottom boundary 
condition is open for in- and outflows. In the simula- 
tions with an initial vertical magnetic field, magnetic 
field lines are forced to be vertical at the bottom bound- 
ary. We also use an alternative configuration of the bot- 
tom boundary condition, that allows the influx of hor- 
izontal magnetic field with a strength of 10 G in the 
convective upflow regions and leaves the outflows unaf- 
fected. This configuration has be en described in more 
detail by Stein & Nordlund (2006) and has been used in 
the simulations of Ste iner et al.l (|2008[ ). 

We use a domain size of 12 x 12 Mm 2 in the horizon- 
tal directions, and 1.4 Mm in the vertical, resolved by 
480 x 480 and 100 grid cells, respectively. The domain 
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Fig. 1. — Vertical component of the magnetic field measured at the photospheric level in four simulations that differ in their initial field 
strength and geometry. Top left: horizontal magnetic field of 10 G advected through the bottom boundary of the numerical domain; top 
right: 200 G initial vertical magnetic field; bottom left: mixed-polarity initial vertical 200 G magnetic field, bottom right: 10 G initial 
vertical magnetic field simulation. The images have been taken approximately 800 s after the magnetic field was introduced into the domain. 



is positioned in such way that the continuum formation 
level is located approximately 800 km above the bottom 
boundary. We use this setup to perform a series of sim- 
ulations with 10 G and 200 G initial uniform vertical 
magnetic fields, as well as a simulation with initial ver- 
tical magnetic field of mixed (positive-negative) polarity 
with a net strength o f 200 G, similar to that used by 
I Cameron et al.l (|2011f ). A simulation with a 10 G hor- 
izontal magnetic field, advected with upflows into the 
domain though its bottom boundary, is also carried out. 

Starting from a non-magnetic photospheric convec- 
tion snapshot, we introduced these magnetic fields and 
boundary conditions in the domain and recorded se- 
quences of magneto-convection snapshots over a few 
granular lifetimes. 

In Fig. [TJ vertical magnetic field distribution B y (x, z) 
in the domain at the photospheric level is shown for four 
simulation runs of different field strength and geometry. 
The top left panel shows that the interaction of hori- 
zontal magnetic field influx with convective flows does 
not form strong intergranular magnetic flux concentra- 
tions, in contrast to the 10 G and 200 G initial verti- 
cal uniform magnetic field (top right and bottom right). 
Small-scale loop-like structures with magnetic fields of 
about 10 — 30 G are formed instead due to the advection 
of magnetic field by granular upflows. The area cov- 
ered by magnetic field concentrations in the 10 G verti- 
cal case is significantly lower compared to that for 200 G, 



and the magnetic field strength in these concentrations is 
also somewhat lower. In the bottom left panel of Fig. [T] 
the distribution of the vertical component of magnetic 
field for the mixed-polarity simulation is shown. Here, 
strong magnetic flux concentrations with strengths of up 
to about 1.2 kG are initially formed. However, the net 
magnetic flux decreases during the simulation as a re- 
sult of oppos ite-polarity magnetic flux reconnection and 
cancellation (jCameron et al. 2011), and due to the mag- 
net ic diffusivity in this relativ ely low-resolution run. (See 
e.g. lVogler fc Schu ssler 2007, for a discussion of magnetic 
field generation and decay speed dependence on the nu- 
merical grid resolution and effective magnetic Reynolds 
number). 

The 10 G initial vertical uniform magnetic field and 
10 G horizontal magnetic field simulations show signifi- 
cant differences in the evolution of their average quanti- 
ties. This is plotted in Fig.[2j where the time evolution of 
the modulus of horizontal components of magnetic field 
is shown for the two cases. In the horizontal magnetic 
field simulation, it takes approximately 15 minutes for 
the magnetic field to propagate from the bottom bound- 
ary of the domain to the photospheric level (marked by 
y = in the left plot of F ig. [2]). In agreement with 
the findings of Steine r~et al.l ((2008), the horizontal mag- 
netic field is accumulated at a height of approximately 
0.15 Mm above the average continuum formation level 
y = 0. This is visible as a small increase in the average 
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Fig. 2. — Time evolution of average horizontal magnetic field strength as a function of height. Left panel: 10 G horizontal magnetic field 
advected through the bottom boundary; right panel: a 10 G initially vertical magnetic field simulation. 




Fig. 3. — Time average of the variation in the mean horizontal 
magnetic field with height. The dashed curve corresponds to the 
10 G initial vertical uniform magnetic field setup, while the black 
solid curve is the 10 G horizontal magnetic field simulation. 



horizontal magnetic field at this height. A similar plot 
for the 10 G vertical field is shown in the right panel of 
Fig. [21 In the latter case, the vertical magnetic field is 
redistributed by the convective flows within the first 5 
minutes of the simulation, and the horizontal component 
of the magnetic field appears. There is also no clear layer 
where the horizontal magnetic field is accumulated. This 
is further demonstrated in Fig. [3j where we plot the time 
averages of the dependence of the horizontal magnetic 
field as a function of height. The figure shows a local 
enhancement of the horizontal magnetic field strength at 
y = 0.15 Mm for 10 G horizontal field simulation, while 
for the vertical magnetic field simulation there is almost 
no enhancement. 

3. POYNTING FLUX ANALYSIS 
We define the MHD Poynting flux vector as 
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B x (v x B) , 



(1) 



where B is the magnetic field vector and v is velocity 
field vector. We are interested in the vertical component 
of the Poynting flux vector F y (note that, in our notation, 
the y vector component corresponds to the vertical, while 
x and z correspond to the horizontal directions) which is 



defined as 



Bl) - B y (v x B x + v z B z )) = F;+F*. 

(2) 

In Eq. ([2]), Fy corresponds to the part of the vertical 
component of the Poynting flux vector produced by ver- 
tical motions, and represents the transport of horizontal 
magnetic field with vertical velocity, while Fy is the part 
of the vertical component, generated by horizontal mag- 
netised plasma motions. 

The F y , Fy and Fy components of the Poynting flux 
vector for 10 G vertical and 10 G horizontal magnetic 
field simulations, measured at 500 km height above the 
average continuum formation level in the simulation do- 
main, are shown in the left, middle and right panels of 
Fig.HJ respectively. For both simulations, the total verti- 
cal (y) component of Poynting flux is generally positive. 
A positive vertical Poynting flux is therefore a general 
property of the simulated photosphere and does not de- 
pend on the initial magnetic field geometry. However, 
our simulations show that the average F y is about 10 
times smaller for the case of horizontal magnetic field 
compared to the 10 G initial vertical magnetic field. This 
finding holds for the latest stages of the simulation which 
have achieved a state of equilibrium. 

The middle panel of Fig. |4] shows the amount of the ver- 
tical Poynting flux, generated by the vertical motions. It 
is evident from the plot that the horizontal magnetic field 
simulation produces an average vertical Poynting flux 
which is positive, while the vertical magnetic field does 
not. This confirms the presence of the mechanism of ver- 
tical Poynting flux production by advection of horizon- 
tal ina^netic_J^din_the upward direction, as described 
bv lSteiner et al.l ([20081 ). However, as shown in the right 
panel of Fig.|4j only about 5% of the total vertical Poynt- 
ing flux is produced by the upward advection of hori- 
zontal magnetic fields. This comparison clearly shows 
that it is the strong vertical magnetic fields in the upper 
photosphere that are mainly responsible for the verti- 
cal Poynting f lux pro duction. As previously shown by 
Shel vag et al.l (|2011b[ h this flux is produced by the inter- 
action of vortex flows in the intergranular lanes with the 
expanding photospheric magnetic concentrations. The 
positive sign (outward direction) of Poynting flux is not 
connected either to the initial magnetic field direction, 
or to the way the magnetic field is introduced into the 
domain. To demonstrate this further, in Fig. [5] we com- 
pare the Poynting flux produced with an initial positive 
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Fig. 4. — The time dependence of the total vertical Poynting flux (left) and its components Fy (middle) and Fy (right) for the 10 G 
vertical (red) and 10 G horizontal simulations. Solid lines in the plots correspond to the time dependences, obtained by smoothing with 
a boxcar window of width 5 minutes, while dotted lines are the non-smoothed dependences. Note that the fluxes for the 10 G horizontal 
magnetic field simulations are multiplied by a factor of 5 to allow their comparison with the 10 G vertical simulations. 
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Fig. 5. — Same as Fig. [2]for 200 G initial vertical and 200 G initial mixed polarity magnetic field simulations. 



vertical 200 G field with that of a mixed-polarity ver- 
tical 200 G field. The dependence of the total vertical 
component of the Poynting flux with time (left panel 
of Fig. [5]) shows that after the first 5 minutes, during 
which the magnetic field is advected into the intergran- 
ular lanes, the behaviour of the vertical magnetic field 
simulation stabilises, while the flux emerging from the 
mixed-polarity simulation begins to decay. In contrast 
to the weaker, 10 G vertical uniform field simulation, the 
vertical Poynting flux generated by vertical motions is 
now negative, suggesting larger amount of stronger hor- 
izontal fields located in the vertical downflow regions. 
During the first 20 minutes of the simulation, the mixed- 
polarity model also shows a negative Poynting flux pro- 
duced by vertical motions (middle panel of Fig. [5]). How- 
ever, after this stage, when a significant amount of verti- 
cal magnetic field has cancelled out, the horizontal small- 
scale magneti c loop-like struct ures in the granular upflow 
regions ( Ca meron et al.l |2011[ ) start producing positive 
Poynting flux, similar to the 10 G horizontal field model. 
As the right panel of Fig. [5] shows, both models generate 
vertical positive Poynting flux due to horizontal motions 
in strong vertical magnetic field concentrations, which is 
about two orders of magnitude larger than the Poynting 
flux produced by vertical motions. 

It has to be noted that the height of 500 km for the 
measurement of the Poynting flux has been chosen ar- 
bitrarily. The results presented remain valid everywhere 
where the Poynting flux is positive, far away from the 
upper boundary and in a wide range of heights starting 
from about 200 km above the photospheric level. Thus, 



the influence of the upper boundary condition is consid- 
ered to be small on these measurements. 

In general, the 10 G simulations show notably more 
turbulent behaviour of Poynting flux compared to the 
200 G models. While the latter shows the oscillations 
in F y , the former demonstrates a more random depen- 
dence of F y with time. This is due to the presence of a 
large relative amount of weak, turbulent granular fields 
compared to the strong magnetic field model, where the 
major part of magnetic flux is concentrated in the inter- 
granular lanes. 

It is also interesting to note that the Poynting flux in 
the vertical direction does not scale simply as the square 
of the magnetic field strength, as may be suggested from 
Eq. (j2j) . While the ratio of squares of the mean magnetic 
field in the 10 and 200 G simulations is 400, the ratio 
of the average vertical Poynting fluxes in these models is 
about 10. This is due to the magnetic field in the inter- 
granular lanes, which does not increase linearly with in- 
creasing initial field. However, it tends to occupy a larger 
area of the intergranular lanes, thus increasing the fill- 
ing factor. This process continues until the intergranular 
lanes are fully filled with the magnetic field. Any further 
increases in the initial magnetic field strength lead to the 
appe arance of a dynamically d ifferent sunspot umbra re- 
gion ([Schussler fc Vo gler 2006|). 

4. CONCLUSIONS 

In this paper we carried out a comparative analysis 
of the Poynting flux generated by simulated magnetised 
solar photospheres with magnetic fields that differ in 
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strength and geometry. We have shown the presence of 
two different types of vertical Poynting flux generation: 
one connected to the advection of horizontal magnetic 
field by upflows in the photosphere, while the other is 
due to horizontal motions of plasma in strong intergran- 
ular magnetic flux concentrations. Using different initial 
magnetic field strength photospheric magneto-convection 
simulations, mixed polarity region simulations and sim- 
ulations with the horizontal magnetic field advected 
through the bottom domain boundary, we showed that 
the latter mechanism is the preferred one, independent 
of both how the magnetic field is introduced in the sim- 
ulated photosphere, and the direction of initial magnetic 
field. We demonstrated clearly that the outward flux of 
electromagnetic energy (Poynting flux) is a general prop- 
erty of the simulated photosphere and does not depend 
on the way how the magnetic field is introduced in the 
domain. 

The initial magnetic field configurations presented in 
this paper can be related to regions of the solar pho- 



tosphere. Our weak magnetic field (10 G) simulations 
correspond to a case of an extremely weakly magnetised 
solar granulation, while the 200 G uniform and opposite- 
polarity models allow production of stronger horizontal 
fields and represent the evolution of unipolar or bipolar 
regions in proximity to sunspots. The results shown sug- 
gest that both magnetically active and non-active granu- 
lation deliver electromagnetic energy to the upper layers 
of solar atmosphere. While simulated active regions with 
strong vertical magnetic fields generate 5 times more 
Poynting flux compared to non- active regions, the area 
covered by non-magnetic granulation can be significantly 
larger than the one covered by active regions, depending 
on the global solar activity level. 
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